Power Law Blinking Quantum Dots: Stochastic and Physical Models 
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We quantify nonergodic and aging behaviors of nanocrystals (or quantum dots) based on stochastic 
model. Ergodicity breaking is characterized based on time average intensity and time average 
correlation function, which remain random even in the limit of long measurement time. We argue 
that certain aspects of nonergodicity can be explained based on a modification of Onsager's diffusion 
model of an ion pair escaping neutralization. We explain how diffusion models generate nonergodic 
behavior, namely a simple mechanism is responsible for the breakdown of the standard assumption of 
statistical mechanics. Data analysis shows that distributions of on and off intervals in the nanocrystal 
blinking are almost identical, ip±{r) cx A±r~ < ' 1+a± ' > with A+ ~ A- and a+ ~ a_ = a and a ~ 0.8. 
The latter exponent indicates that a simple diffusion model with a = 0.5 neglecting the electron-hole 
Coulomb interaction and/or tunneling, is not sufficient. 



INTRODUCTION 



Single quantum dots when interacting with a continu- 
ous wave laser field blink: at random times the dot turns 
from a state on, in which many photons are emitted to a 
state off in which no photons are emitted. While stochas- 
tic intensity trails are found today in a vast number of 
single molecule experiments, the dots exhibit statistical 
behavior which seems unique. In particular, the dots 
exhibit power law statistics, aging, and ergodicity break- 
ing. While our understanding of the Physical origin of 
the blinking behavior of the dots is not complete, several 
physical pictures have emerged in recent years, which 
explain the blinking in terms of simple Physics. Here 
we will review a diffusion model which might explain 
some of the observations made so far. Then we analyze 
the stochastic properties of the dots, using a stochas- 
tic approach. In particular we review the behaviors of 
the time and ensemble average intensity correlation func- 
tions. Usually it is assumed that these two objects are 
identical in the limit of long times, however this is not 
the case for the dots. 



II. PHYSICAL MODELS 

A typical fluorescence intensity trace of a CdSe quan- 
tum dot, or nanocrystal (NC), overcoated with ZnS (in 
short, CdSe-ZnS NC) under continuous laser illumina- 
tion is shown in Fig. ^ From this Figure we learn, that 
roughly, the intensity jumps between two states - on and 
off. Some of the deviations from this digital behavior can 
be attributed to fluctuating non-radiative decay channels 
due to coupling to the environment, and also to time bin- 
ning procedure 0.Q [H| an d see also Q. Data analysis of 
such time trace is many times based on distribution of on 
and off times. Defining a threshold above which the NC 
is considered in state on and under which it is in state off, 
one can extract the probability density functions V±( T ) 
of on and ip-(T) of off times. Surprisingly these show a 
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Figure 1: Intensity fluctuations in a CdSe-ZnS NC under con- 
tinuous laser illumination at room temperature. Dotted hor- 
izontal line was selected as a threshold to divide off and on 
states. 



power-law decay ip±{r) cx r~ 1 ~ a± , as shown in Fig. |2 
A summary of different experimental exponents is pre- 
sented in Table HI indicating such a power-law decay in 
most cases. In some cases a + « a_ and the exponents 
are close to 1/2. In particular, Brokmann et al. @ mea- 
sured 215 CdSe-ZnS NCs and found that all are statisti- 
cally identical with a+ = 0.58 ±0.17, a_ w 0.48 ±0.15 so 
that a+ w a_ w 0.5. Note that most of the uncertainty 
in the values of the exponents can be attributed simply to 
statistical limitations of data analysis (see also Sec- 
tion El below). Shimizu et al. found that in the limit 
of low temperature and weak laser fields a+ « a_ w 0.5. 
The fact that in many cases a± < 1, leads to interest- 
ing statistical behavior, for example ergodicity breaking, 
and aging. We will discuss these behaviors in Sec. III. A 
physical model for blinking was suggested by Efros and 
Rosen Briefly the on and off periods correspond to 
neutral and charged NCs respectively. Thus the on/off 
trace teaches us something on elementary charging mech- 
anism of the dot. The difficulty is to explain the power 
law distributions of on and off times, or in other words 
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Group 


Material 


No. 


Radii, nm 


Temp., K 


Laser Intensity, 


a+ 


OL- 


Verberk et al. [4] 


CdS 


1 


2.5 


1.2 




e~ at 


0.65(0.2) 


Brokmann et al. [5J 


CdSe-ZnS 


215 




300 




0.58(0.17) 


0.48(0.15) 


Shimizu et al. [6] 


CdSe-ZnS, CdSe, CdTe 


>200 


1.5, 2.5 


300, 10 


0.1-0.7 


0.5(0.1), cutoff 


0.5(0.1) 


Kuno et al. [JJ 


CdSe-ZnS 


~ 200 


1.7-2.9 


300-394 


0.24-2.4 




0.5-0.75 


Kuno et al. [8J 


CdSe-ZnS 


>300 


1.7-2.7 


300 


0.1-100 


0.9(0.05) 


0.54(0.03) 


Kuno et al. [9J 


InP 


~ 30 


1.5 


300 


0.24 


1.0(0.2) 


0.5(0.1) 


Cichos et al. flO] 


Si 








1.8, 6.5 


1.2(0.1) 


0.3, 0.7 


Hohng and Ha [11] 


CdSe-ZnS 


~ 1000 










0.94-1.10 


Miiller et al. [12] 


CdSe-ZnS 




4.4 (core) 


300 


0.025 


0.55 


0.05, 0.25 


van Sark et al. [13] 


CdSe-ZnS 


41 


~ 3.7 


300 


20 


~ 1.2, ~ 0.7 


~ 0.2, ~ 0.4 


Kobitski et al. [14] 


CdSe 




3.6 




0.04-0.38 


0.97-0.66 


0.42-0.64 



Table I: Summary of experimental exponents for on (a+) and off (a_) time distributions for various single NCs under different 
experimental conditions. Notice that Verberk et al. use uncapped NCs, while other measurement consider capped NCs, hence 
exponential distribution on times is found only for uncapped dots. Hohng and Ha used CdSe-ZnS NCs coated with streptavidin 
which might alter the exponent a~. 
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Figure 2: Distributions i/>±(t) of on and off times for the NC, 
whose intensity trajectory is shown in Fig. U The straight 
line is the fit to the off time distribution. 



why should the time the charge occupies the NCs follow 
power law behavior? 

Two types of models were suggested, a diffusion ap- 
proach and a random trap model. The measurements 
of Dahan's and Bawendi's groups 0,0, which show the 
universal power law a± = 0.5, are consistent with the 
diffusion model (see details below). The fact that all 
dots are found to be similar jalseem not consistent with 
models of quenched disorder 0, 0, O] since these sup- 
port the idea of a distribution of a±. However, some 
experiments show deviations from the a+ « a_ rj 0.5 
and might support the distribution of a± . It is possible 
that preparation methods and environments lead to dif- 
ferent mechanisms of power law blinking, and different 
exponents ^^|. More experimental work in this direc- 



tion is needed, in particular, experimentalists still have 
to investigate the distribution of a± , and show whether 
and under what conditions are all the dots statistically 
identical. Below we discuss the diffusion model; differ- 
ent aspects of the tunneling and trapping model can be 

found in HE E3. 

As discussed at length by Shimizu et al. 0], the on 

time distributions show temperature and laser power de- 
pendencies, e.g. exponential cutoffs of power law behav- 
ior. Although no direct observations of cutoffs in the off 
time distribution was reported, ensemble measurements 
by Chung and Bawendi [j] demonstrate that there should 
be such a cutoff as well, but at times of the order of tens 
of minutes to hours. Our analysis here, employing the 
power law decaying distributions, is of course applicable 
in time windows where power law statistics holds. 



A. Diffusion model 

We note that the simplest diffusion controlled chem- 
ical reaction A + B ^ AB, where A is fixed in space, 
can be used to explain some of the observed behavior 
on the uncapped NCs. As shown by the group of Orrit 
[1] such dots exhibit exponential distribution of on times 
and power law distribution of off times. The on times 
follow standard exponential kinetics corresponding to an 
ionization of a neutral NC (denoted as AB). A model 
for this exponential behavior was given already in flfl ]. 
Clearly the experiments of the group of Orrit, show that 
the capping plays an important part in the the blink- 
ing, since capped NCs exhibit power law behavior both 
for the on and off times. We will return to capped dots 
later. 

Once the uncapped NC is ionized (A + B state) we 
assume the ejected charge carrier exhibits a random walk 
on the surface of the NC or in the bulk. This part of 
the problem is similar to Onsager's classical problem of 
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an ion pair escaping neutralization (see e.g., [Isl 
The survival probability in the off state for time t, S- (t) 
is related to the off time distribution via S-(t) = 1 — 
f*t/>-(T)dT, or 



t/>-(t) = - 



dS-(t) 
dt ' 



(1) 



It is well known that in three dimensions survival prob- 
ability decays like t -1 / 2 , the exponent 1/2 is close to the 
exponent often measured in the experiments. In infinite 
domain the decay is not to zero, but the 1/2 appears in 
many situations, for finite and infinite systems, in com- 
pletely and partially diffusion controlled recombination, 
in different dimensions, and can govern the leading be- 
havior of the survival probability for orders of magnitude 
in time [lflll^ l2l| . In this picture the exponent 1/2 does 
not depend on temperature, similar to what is observed 
in experiment. We note that it is possible that instead of 
the charge carrier executing the random walk, diffusing 
lattice defects which serve as a trap for charge carrier are 
responsible for the blinking behavior of the NCs. 

A long time ago, Hong, Noolandi and Street jjjjjj inves- 
tigated geminate electron-hole recombination in amor- 
phous semiconductors. In their model they included 
the effects of tunneling, Coulomb interaction, and dif- 
fusion. Combination of tunneling and diffusion leads to 
a S(t) oc t -1 / 2 behavior. However, when the Coulomb 
interactions are included in the theory, deviations from 
the universal t -1 ' 2 law, are observed. For example in the 
the analysis of photoluminescence decay in amorphous 
Si:H, as a function of temperature. 

Coulomb interaction between the charged NCs and the 
ejected electron seems to be an important factor in the 
Physics of NCs. The Onsager radius is a measure of the 
strength of the interaction 



ror, 



k h Te 



(2) 



Krauss and Brus |23] measured the dielectric constant of 
CdSe dots, and found the value of 8. Hence, at room 
temperature we find rons — 70A (however, note that the 
dielectric constant of the matrix is not identical to that of 
the dot). Since the length scale of the dots is of the order 
of a few nanometers, the Coulomb interaction seems an 
important ingredient of the problem. This according to 
the theory in j^j is an indication of possible deviations 
from the universal 1/2 power law behavior. It is also an 
indication that an ejected electron is likely to return to 
the dot and not escape to the bulk (since the force is 
attractive). In contrast, if the Onsager radius is small, 
an ejected electron would most likely escape to the bulk, 
leaving the dot in state off forever (i.e. Polya theorem 
in three dimensions). Unfortunately, currently there is 
not sufficient experimental data to determine in more 
qualitative ways if, Onsager type of model can be used 
to explain the observed data. As in standard geminate 
recombination processes, the dependence of blinking on 
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Figure 3: On and off states for NCs, following Q. 



temperature, dielectric constant of the dot and of the 
matrix and on external driving field, might yield 
more microscopical information on the precise physical 
mechanism of the fascinating blinking behavior. 

One of the possible physical pictures explaining blink- 
ing of capped NCs can be based on diffusion process, 
using a variation of a three state model of Verberk et al. 

As mentioned above, for this case power law distri- 
bution of on and off times are observed. In particular, 
neutral capped NC will correspond to state on (as for 
uncapped NCs). However, capped NC can remain on 
even in the ionized state - see Fig. Verberk et al. as- 
sume that the ionized capped NC can be found in two 
states: (i) the charge remaining in the NC can be found 
in center of NC (possibly a de-localized state), (ii) charge 
remaining in the NC can be trapped in vicinity of cap- 
ping. For case (i) the NC will be in state off, for case (ii) 
the NC will be in state on. Depending on exact location 
of this charge, the fluorescence intensity can vary. The 
main idea is that the rate of Auger nonradiative recom- 
bination ^3] of consecutively formed electron-hole pairs 
will drop for case (ii) but not for case (i). We note that 
capping may increase effective radius of the NC, or pro- 
vide trapping sites for the hole (e.g., recent studies by 
Lifshitz et al. [131 demonstrate that coating of NCs cre- 
ates trapping sites in the interface). Thus the off times 
occur when the NC is ionized and the hole is close to the 
center, these off times are slaved to the diffusion of the 
electron. While on times occur for both a neutral NC 
and for charged NC with the charge in vicinity of cap- 
ping, the latter on times are slaved to the diffusion of 
the electron. In the case of power law off time statistics 
this model predicts same power law exponent for the on 
times, because both of them are governed by the return 
time of the ejected electron. 

Beyond nanocrystals, we note that fluorescence of sin- 
gle molecul es [2^1 and of nanoparticles diffusing through a 
laser focus |26| , switching on and off of vibrational modes 
of a molecule [13|, openin g-cl osing behavior of certain 
single ion channels [20L Hh. I2i3| . motion of bacteria 
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Figure 4: Schematic temporal evolution of the dichotomous 
intensity process. 



deterministic diffusion in chaotic systems 13 U. the sign of 
magnetization of spin systems at criticality l32|. and oth- 
ers exhibit power law intermittency behavior [33]. More 
generally the time trace of the NCs is similar to the well 
known Levy walk model piij . Hence the stochastic theory 
which we consider in the following section is very general. 
In particular we do not restrict our attention to the ex- 
ponent 1/2, as there are indications for other values of a 
between and 1, and the analysis hardly changes. 



III. STOCHASTIC MODEL AND DEFINITIONS 

The random process considered in this manuscript, is 
shown in Fig. 31 The intensity I(t) jumps between two 
states I(t) = +1 and I(t) = 0. At start of the measure- 
ment t = the NC is in state on: 1(0) = 1. The sojourn 
time n is an off time if i is even, it is an on time if i 
is odd (see Fig. BJ. The times t% for odd [even] i, are 
drawn at random from the probability density function 
(PDF) ip+(t), [ip-(t)], respectively. These sojourn times 
are mutually independent, identically distributed random 
variables. Times are cumulative times from the pro- 
cess starting point at time zero till the end of the i'th 
transition. Time T' on Fig. 0] is the time of observation. 
We denote the Laplace transform of ip± (t) using 



= / i>±{t)e- st dt. 



(3) 



In what follows we will investigate statistical properties 
of this seemingly simple stochastic process. In particular 
we will investigate the correlation function of this pro- 
cess. In experiment correlation functions are used many 
times to characterize intensity trajectories. The main ad- 
vantage of the analysis of correlation functions, if com- 
pared with PDFs of on and off times, is that in former 
case there is no need to introduce the intensity cutoff. 
Correlations functions are more general than on and off 
time distributions. Besides, correlation functions exhibit 
aging, and ergodicity breaking, which are in our opinion 
interesting. 

We will consider several classes of on/off PDFs, and 
classify generic behaviors based on the small s expansion 
of ip± (s). We will consider: 



(i) Case 1 PDFs with finite mean on and off times, 
whose Laplace transform in the limit s — > satisfies: 



1p±(s) = 1 — ST± 



(4) 



Here t + (t_) is the average on (off) time. For example 
exponentially distributed on and off times, 



i 



1 + ST± 



(5) 



belong to this class of PDFs. 

(ii) Case 2 PDFs with infinite mean on and off times, 
namely PDFs with power law behavior satisfying 



%p± oc t 



-1 — <X-\ 



ot- < a + < 1, 



(6) 



in the limit of long times. The small s behavior of these 
family of functions satisfies 



^±(s) = 1 - A±s a± 



(7) 



where A± are parameters which have units of time" . We 
will also consider cases where on times have finite mean 
(a + = 1) while the off mean time diverges (a- < 1) 
since this situation describes behavior of uncapped NC 
(see also 0). 

(iii) Case 3 PDFs with infinite mean with a + = a_ = a 



i>±(s) = 1 - A±s a 



(«) 



As mentioned Brokmann et al. @ report that for CdSe 
dots, a+ = 0.58 ± 0.17, and a- = 0.48 ± 0.15, hence 
within error of measurement, a ~ 0.5. 

Standard theories of data analysis, usually use the er- 
godic hypothesis and a time average of a process is re- 
placed with an average over an ensemble. The simplest 
time average in our case is the time average intensity 



7 if jm 



(9) 



In the limit of long times and if ergodic assumption holds 
I = {I), where (7) is the ensemble average. As usual we 
may generate many intensity trajectories one at a time, 
to obtain ensemble averaged correlation function 



C(t,t') = {I(t)I(t + t% 



(10) 



and the normalized ensemble averaged correlation func- 
tion 

9 { > ] ~ (/(*)> </(*+*)) (/(*)) <j(t +*)>' 1 ] 

From a single trajectory of I(t), recorded in a time in- 
terval (0, T'), we may construct the time average (TA) 
correlation function 
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In single molecule experiments, the time averaged corre- 
lation function is considered, not the ensemble average. 
However, it is many times assumed that the ensemble 
average and the time average correlation functions are 
identical. For nonergodic processes Cta(T' ' ,t r ) ^ C(t,t') 
even in the limit of large t and T 1 , Moreover for noner- 
godic processes, even in the limit of T" — > oo, Cta(T', t') 
is a random function which varies from one sample of 
I(t) to another. The ensemble-averaged function C(t,t') 
of the considered process is non-stationary, i.e., it keeps 
its dependence on t even when t — > oo. This is known as 
aging. It follows then from Eq. iJEJ that (C TA (T', t')) = 

C{t,t>)dt/{T -t')^C{t,t<). 



IV. AGING 

Consider the ensemble averaged correlation function 
C(t,t') — (I(t + t')I(t)}, For processes with finite micro- 
scopical time scale, which exhibit stationary behavior, 
one has C(t,t') — f(t'). Namely the correlation func- 
tion does not depend on the observation time t. Aging 
means that C(t, t') depends on both t and t' even in the 
limit when both are large (IE H3 ■ Simple aging behavior 
means that at the scaling limit C(t, t') = f(t'/t), which is 
indeed the scaling in our Case 3; in Case 2 below we find 
such a scaling for g^(t, t'), while C(t, t') will scale differ- 
ently. Aging and non-ergodicity are related. In our mod- 
els, when single particle trajectories turn non-ergodic, the 
ensemble average exhibit aging. Both behaviors are re- 
lated to the fact that there is no characteristic time scale 
for the underlying process. 



A. Mean Intensity of on-off process 

The ensemble averaged intensity (I(t)) for the process 
switching between 1 and and starting at 1 is now con- 
sidered, which will be used later. In Laplace t — > s space 
it is easy to show that 



I( 8 ) 



1 



(13) 



The Laplace s — > t inversion of Eq. l(T3l) yields the mean 
intensity {I(t)). Using small s expansions of Eq. lfl3ll . 
we find in the limit of long times 



(A*)> 



T++T- 

A + t a -- a + 
A_r(l+a_-a + ) 

A± 

A++A- 



case 1 



case 2 



case 3. 



(14) 



If the on times are exponential, as in Eq. (0 then 



1(8) 



1 + ST -i 



1>-(*Y 



(15) 



This case corresponds to the behavior of the uncapped 
NCs. The expression in Eq. ^J, and more generally, 
the case a- < a + = 1 leads for long time t to 



r+t c 



A_T(a_) 



(16) 



For exponential on and off time distributions Eq. (0, 
we obtain the exact solution 



(Ht)) 



r_ exp 



-t 



(17) 



The average intensity does not yield direct evidence 
for aging, because it depends only on one time variable, 
and one has to consider a correlation function to explore 
aging in its usual meaning. 

Remark For the case a + < a_ < 1, corresponding to 
a situation where on times are in statistical sense much 
longer then off times, (I(t)) ~ 1. 



B. Aging Correlation Function of on-off process 

The ensemble averaged correlation function C(t,t') — 
(I(t)I(t + t')) was calculated in j^j. Contributions to 
the correlation function arise only from trajectories with 
I(t) = 1 and I(t + t') = l, yielding 



C(t,u) 



f t (u = 0,+)-f t (u,+) 



4>-(u) l-$ + (u) 



+ ft(u,+y 

u 1 — 4>-(u)ip+(u) 

where u is the Laplace conjugate of t' and 

$+(s) - tp+(u) 



(18) 



/.(«,+) 



(u-s) 1 - 4> + (s)ip-(s) 



(19) 



where s is the Laplace conjugate of t. We note that 
f s (u, +) is the double Laplace transform of the PDF of 
the so called forward recurrence time. This means that 
after the aging of the process in time interval t, the statis- 
tics of first jump event after time t will generally depend 
on the age t. However, a process is said to exhibit aging, 
only if the statistics of this first jump depend on t even 
when this age is long. In particular if the microscopical 
time scale of the problem is infinite, no matter how big 
is t the correlation function still depends on the age (see 
details below). The first term in Eq. fl8jl is due to tra- 
jectories which were in state on at time t and did not 
make any transitions (i.e. the concept of persistence), 
while the second term includes all the contributions from 
the trajectories being in state on at time t and making 
an even number of transitions [331 . 
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C. Case 1 

For case 1 with finite r + and t_, and in the limit of 
long times t, we find 

lim C(t, u) = 

t — >oo 



1 — lp+i u ) 1 — '0-( it ) 



u r+ 



T+1t 



1 - ^(u)tp+(u) 



(20) 



This result was obtained by Verberk and Orrit 39] and it 
is seen that the correlation function depends asymptot- 
ically only on t' (since u is Laplace pair of t'). Namely, 
when average on and off times are finite the system does 
not exhibit aging. If both ^p+(t) and 4>-(t) are exponen- 
tial then the exact result is 



C(t,t>) = 



T- exp 



-t 



T- exp 



- T + 



and C(t,t') becomes independent of t exponentially fast 
as t grows. 



D. Case 2 

We consider case 2, however limit our discussion to 
the case a + — 1 and a_ < 1. As mentioned this case 
corresponds to uncapped NCs where on times are expo- 
nentially distributed, while off times are described by 
power law statistics. Using the exact solution Eq. l|T8jl 
we find asymptotically, when both t and t' are large: 



(21) 



Unlike case 1 the correlation function approaches zero 
when t — > oo, since when t is large we expect to find 
the process in state off. Using Eq. iflfill . the asymptotic 
behavior of the normalized correlation function Eq. 11 111 
is 
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(2) 



(t,f) 



1 



(22) 



We see that the correlation functions Eqs. l(2"Tl l22|l ex- 
hibit aging, since they depend on the age of the process 
t. 

Considering the asymptotic behavior of C(t, t') for 
large t, 

C(t,u) w 



1 r+ 
uA_T(a-)t 1 - 



1 - 



1 - V- 



T + U 



1 



-0_(u)-0+(w) 

(23) 

This equation is similar to Eq. ffity. especially if we 
notice that the "effective mean" time of state off until 
total time t scales as A_i 1-Q -. 

For the special case, where on times are exponentially 
distributed, the correlation function C is a product of 
two identical expressions for all t and t'; 



C{s,u) 



(24) 



where s (u) is the Laplace conjugate of t (t') respectively. 
Comparing to Eq. 11511 we obtain 



C(t,t') = (I(t))(I(t')}, 
and for the normalized correlation function 



5 (2) (M') = 



(i(t + t')Y 



(25) 



(26) 



Eqs. 11261 12511 are important since they show that mea- 
surement of mean intensity (1(f)) yields the correlation 
functions, for this case. While our derivation of Eqs. il26l 
l25ll is based on the assumption of exponential on times, 
it is valid more generally for any ip+(t) with finite mo- 
ments, in the asymptotic limit of large t and t' . To see 
this note that Eqs. plllTfi| yield C(t,t') - (I(t)) (1(f)). 

In Fig. we compare the asymptotic result (12H with 
exact numerical double Laplace inversion of the corre- 
lation function. We use exponential PDF of on times 
ip+{s) = 1/(1 + s), and power law distributed off times: 
'ijj-(s) = V>-(s) = 1/(1 + s QA ) corresponding to a- = 0.4. 
Convergence to asymptotic behavior is observed. 

Remark For fixed t the correlation function in Eq. 
plf exhibits a (£')"" _1 decay. A (t') 01 -- 1 decay of an in- 
tensity correlation function was reported in experiments 
of Orrit's group for uncapped NCs (for that case 
q_ = 0.65 ± 0.2). However, the measured correlation 
function is a time averaged correlation function Eq. l(T2|l 
obtained from a single trajectory. In that case the cor- 
relation function is independent of t, and hence no com- 
parison between theory and experiment can be made yet. 



E. Case 3 



We now consider case 3, and find 



C(t,t')=P + -P +P J-^B 



l + t/t 



7 ;1- a,a ) , 



(27) 



where 



P± = 



At 
A+ + A- 



Figure 5: Exact C(t,t') for Case 2: exponential on times 
and power law off times with a_ = 0.4. We use ip+(s) = 
1/(1 + s) and $-(s) = 1/(1 + s 0,4 ) and numerically obtain 
the correlation function. For each curve in the figure we fix 
the time t. The process starts in the state on. Thick dashed 
straight line shows the asymptotic behavior Eq. 12111 . For 
short times (t' < 1 for our example) we observe the behavior 
C(t,t') ~ C(t,0) = (I(t)), the correlation function is flat. 



following from Eq. l|Tlj) . and where 

B(z;a,b)= [ z a-1 (l - xf^dx 
Jo 

is the incomplete beta function. The behavior in this 
limit does not depend on the detailed shape of the PDFs 
of the on and off times, besides the parameters A+/A- 
and a. We note that both terms of Eq. iflSjl contribute to 
Eq. I|27|) . The appearance of the incomplete beta func- 
tion in Eq. I|27|) is related to the concept of persistence. 
The probability of not switching from state on to state 
off in a time interval (t,t + t'), assuming the process is in 
state on at time t, is called the persistence probability. 
In the scaling limit this probability is 



Figure 6: Exact C(t, t') for case 3, when both on and off times 
are power law distributed with a = 0.4. We use ip±(s) = 
1/(1 + s ' 4 ) for different times t increasing from the topmost 
to the lowermost curves. The dots on the left and on the right 
show C(t,0) = (/(*)) and C(t, oo) = (/(*)) /2 respectively. 
The process starts in the state on. 



tails the two limiting behaviors are: 

f P+ £ « 1 

(29) 

Using Eq. 11 411 the normalized intensity correlation func- 
tion is g {2 \t,t') ~ C{t,t')/(P+) 2 . 

In Fig. we compare the asymptotic result 112 7|l with 
exact numerical double Laplace inversion of the correla- 
tion function for PDFs iJj + (s) = -0_(s) = 1/(1 + s 0A ). 
Convergence to Eq. l|27Tl is seen. 

Remark For small t'/t we get flat correlation func- 
tions. Flat correlation functions were observed by Da- 
lian's group for capped NCs. However, the mea- 
sured correlation function is a single trajectory correla- 
tion function Eq. l(T2|l , and hence no comparison between 
theory and experiment can be made yet. 



„ , / \ sin7ra „ / 1 \ 
P (t,t + t')~l —B f — —f ;l-a,a) . (28) 

The persistence implies that long time intervals in which 
the process does not jump between states on and off, 
control the asymptotic behavior of the correlation func- 
tion. The factor P+, which is controlled by the ampli- 
tude ratio A+/A-, determines the expected short and 
long time t' behaviors of the correlation function, namely 
C(oo,0) = Umt^ 00 (7(t)J(t + Q)) = P+ and C(oo,oo) = 
lim t ^ 00 (/(t)/(i + oo)) = {P+f. With slightly more de- 



V. NON ERGODICITY 

Non-ergodicity of blinking quantum dots was first 
pointed out in the experiments of the group of Dahan 
|4Tl |. We begin the discussion of nonergodicity in blink- 
ing NCs by plotting 100 time averaged correlation func- 
tions from 100 NCs in Fig. Clearly, correlation func- 
tions obtained are different. The simplest explanation 
would be that the NCs have different statistical prop- 
erties. However, similar variability is also observed for 
a given NC, when we calculate correlation functions for 
different T" (e.g., H3]). To further illustrate this point, 
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= f / T 



Figure 7: 100 experimental time averaged correlation func- 
tions (one of which is obtained from the signal shown in Fig. 

after "renormalizing" the average on and off intensities to 
be 1 and 0, respectively. Note logarithmic abscissa. 




r = f / T 



Figure 8: Ten typical simulated realizations of Cta for a = 
0.8. 



we generate on a computer the two state process, with 
power law waiting time of on and off times following 
V'-(t) = "0+( r ) = (ff" 1- " for t > 1 (and zero other- 
wise). For each trajectory we calculate its own time av- 
erage correlation. As we show in Fig. [HI the trajectories 
exhibit ergodicity breaking. The most striking feature 
of the figure is that even though trajectories are statis- 
tically identical, the correlation function of the process 
is random, similar to the experimental observation. In 
complete contrast, if we consider a two state process with 
on and off times following exponential statistics, then all 
the correlation functions would be identical, and all of 
them would follow the same master curve: the ensemble 
average correlation function. 

In this section we consider the non-ergodic properties 
of the blinking NCs using a stochastic approach. We as- 
sume that all the NCs are statistically identical in agree- 
ment with 0|, and restrict ourselves to the Case 3. For 
the sake of simplicity we only consider the case when 
distribution of on times is identical to distribution of off 



times, namely a_ = a + = a and A + = A_ = A. Gener- 
alization to A + ^ A- is straightforward [ij. The noner- 
godicity is found only for a < 1, when the mean transi- 
tion time is infinite, and should therefore disappear when 
exponential cutoffs of off and on times become relevant 
0, i-e., when the mean transition times become of the 
order, or less than the experimental time. The described 
model, however, is valid in a wide time window spanning 
many orders of magnitude for the NCs, and is relevant 
to other systems, as mentioned in Section ITT1 



A. Distribution of time averaged intensity 

As mentioned in the introduction, the blinking NCs 
exhibit a non-ergodic behavior. In particular the ensem- 
ble average intensity (I) is not equal to the time average 
I. Of course in the ergodic phase, namely when both the 
mean on and off times are finite, we have (7) = I, in the 
limit of long measurement time. More generally we may 
think about I as a random function of time, which will 
vary from one measurement to another. In the ergodic 
phase, and in the asymptotic limit the distribution of I 
approaches a delta function 



P(I)^S(I -(/}). 



(30) 



The theory of non-ergodic processes deals with the ques- 
tion what is the distribution of P{T) in the non-ergodic 
phase. For the two state stochastic model 



(31) 



where T + is the total time spent in state on. 

A well known example of similar ergodicity breaking 
is regular diffusion, or a binomial random walk on a line. 
The walker starts at the origin and can go left or right 
randomly, at each step. Let the measurement time be 
t, and the position of the random walker be x(t). The 
total time the walker remains on the right of the origin 
x(t) > is T + . The PDF of return time (or of number of 
steps) t to the origin decays as t~ 3 / 2 for large t, so that 
a = 1/2. Two half- axes at both sides of the origin can be 
thought of as the two states, on and off, of the random 
walker. The well-established result is that the fraction I 
of total time spent by the walker on either side, in the 
long time limit is given by the arcsine law (22, Ui| 



PI) 



1 



A main feature of this PDF is its divergence at I — 0, 1, 
indicating that the random walker will most probably 
spend most of its time on one side (either left or right) 
of the origin. In particular the naive expectation that 
the particle will spend half of its time on the right and 
half on the left, in the limit of long measurement time, 
is wrong. In fact the minimum of the arcsine PDF is 
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0.2 0.4 0.6 0.8 1 

T*/T 



Figure 9: The probability density function of I = T + /T for 
the case ip+(i) = ip-(t) <x For the ergodic phase 

a > 1 P(I) is a delta function on (I) = 1/2. In the non- 
ergodic phase 7 is a random function, for small values of a 
the P(I) is peaked on I = and 1 = 1, indicating a trajectory 
which is in state off or on for a period which is of the order 
of measurement time T. 




0.2 0.4 0.6 0.8 1 0.2 0.4 0.6 0.8 1 



C ta' C ta' 

Figure 10: PDF of Cta(T' ,t') for different r = t'/T' and 
a = 0.3. Abscissas are possible values of Cta(T' ,t'). Dia- 
monds are numerical simulations. Curves are analytical re- 
sults without fitting: for r = Eq. 13211 is used (full line), for 
r = 0.01 and 0.1 Eq. I|35fl is used (dashed) and for r = 0.5, 
0.9 and 0.99 Eq. QS) is used (full). 



on I = (I) = 1/2. In other words the ensemble average 
(I) = 1/2 is the least likely event. This result might seem 
counter intuitive at first, but it is due to the fact that 
the mean time for return to the origin is infinite. This 
in turn means that the particle gets randomly stuck on 
x < or on x > for a period which is of the order of the 
measurement time, no matter how long this measurement 
time is. 

In the more general case < a < 1 the distribution of 
/ can be calculated based on the work of Lamperti 
(see also (23), and one finds 

= sin fra) 1^ {l -1)^ 

* 7 2Q + (i-7) 2 " + 2cos(™)7 Q (l-7) a 

(32) 

which is shown in Fig. El When a — > the PDF of I is 
peaked around 1 = and I = 1, corresponding to blink- 
ing trajectories which for most of the observation time T 
are in state off or state on respectively. When a — > 1 , we 
see that l a (I) attains a maximum when I = {I) = 1/2, 
indeed in the ergodic phase a > 1 we obtain as expected a 
delta peak centered on I = 1/2, as we mentioned. There 
exists a critical a c = 0.594611... above (under) which 
l a (I) has a maximum (minimum) on I = 1/2. Note that 
the Lamperti PDF in Eq. i!32|) is not sensitive to the pre- 
cise shapes of the on and off time distributions (besides 
a of course). For situations in which A- ^ A + the sym- 
metry of the Lamperti PDF will not hold. Note that line 
shapes with structures similar to those in Fig. were 
obtained by Jung et al. [Hj in a related problem. Simi- 
lar expressions are also used in stochastic models of spin 
dynamics [45j, and in general, the problem of occupation 




0.2 0.4 0.6 0.8 1 0.2 0.4 0.6 0.8 1 



C ta' C ,a 

Figure 11: PDF of Cta(T' ,t') for different r = t'/T' and 
a = 0.5. Diamonds are numerical simulations. Curves are 
analytical results without fitting: for r = Eq. Q32H is used 
(full line), for r = 0.01 and 0.1 Eq. 13511 is used (dashed) and 
for r = 0.5, 0.9 and 0.99 Eq. Q29 is used (full). 



times, and a related persistence concept, are of a wide 
interest in different fields (ill EfiL l47j . 

Next we extend our understanding of the distribution 
of time averaged intensity to the time averaged correla- 
tion functions defined in Eq. 11211 . 
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the time average intensity between time a and time b > a, 
and 



T = T' — t', 



V 



' 0.2 0.4 0.6 0.8 1 0.2 0.4 0.6 0.8 1 

Figure 12: PDF of Cta(T' ,t') for different r = t'/T' and 
a = 0.8. Diamonds are numerical simulations. Curves are 
analytical results without fitting: for r = Eq. %\2l is used 
(full line), for r = 0.01, 0.1 and 0.5 Eq. is used (dashed) 
and for r = 0.9 and 0.99 Eq. 13St is used (full). If compared 
with the cases a — 0.3 and 0.5, the distribution function 
exhibits a weaker non-ergodic behavior, namely for r = the 
distribution function peaks on the ensemble average value of 
1/2. 



B. Distribution of time averaged correlation 
function 

We first consider the non-ergodic properties of the cor- 
relation function for the case t' = 0. It is useful to define 

2[«,6] - / I(t)dt/(b-a), (33) 

J a 

i 



Using Eq. ((12)) and for t' — the time averaged correla- 
tion function is identical to the time average intensity 



C T a(T,0)=X [0 , t] = (34) 



and its PDF is given by Eq. Figs. EH Eland US for 

the case r = 0, show these distributions for a = 0.3, 0.5 
and 0.8, respectively, together with the numerical results. 

An analytical approach to estimate the distributions 
Pn^ „ it' t.n (z) o f Cta (T\ t') = z for nonzero t' was devel- 
oped in jj33, io] • To treat the problem a non-ergodic mean 
field approximation was used, in which various time av- 
erages were replaced by the time average intensity 2T[o,r]> 
specific for a given realization. For short t' <g; T" the 
result is 



C TA (T',t') 



T 2 

[0,T] 



1-Xr 



[0,T], 



/ \ 1~ a 

\(l-r)I [0T] ) \ na ^ 1 ) 



(l-r)I[ 0T j 



t' < T+ 
t' > T+. 



(35) 



Eq. i'>")li yields the correlation function, however unlike 
standard ergodic theories the correlation function here is 
a random function since it depends on X[o,t]- The distri- 
bution of Cta(T' ',t') is now easy to find using the chain 
rule, and Eqs. (3331 EH- In Figs. UHlHIlandHlwe plot 
the PDF of Cta(T', t') (dashed curves) together with nu- 
merical simulations (diamonds) and find excellent agree- 
ment between theory and simulation, for the cases where 
our approximations are expected to hold r < 1/2. We 
observe that unlike the r = case the PDF of the corre- 
lation function exhibit a non-symmetrical shape. To un- 
derstand this note that trajectories with short but finite 
total time in state on (T + ^ T) will have finite correla- 



tion functions when t' — 0. However when t 1 is increased 
the corresponding correlation functions will typically de- 
cay very fast to zero. On the other hand, correlation 
functions of trajectories with T + ~ T don't change much 
when t' is increased (as long as t' ^ T + ). This leads to 
the gradual nonuniform shift to the left, and "absorption" 
into Cta(T' ,t') — 0, of the Lamperti distribution shape, 
and hence to non-symmetrical shape of the PDFs of the 
correlation function whenever r^O. 

We now turn to the case T <Z^t' . Then 



Cta{T', t 1 ) ~ l[o t T]I[t',T'] ■ (36) 
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In the limit t'/T' — > 1 this yields 

Pc TA (T',t')(z)~[£ a (z) + 6(z)}/2, (37) 

which is easily understood if one realizes that in this limit 
2u',t'] m Eq. i|36|) is either or 1 with probabilities 1/2, 



and that the PDF of T [0 ,T] is Lamperti's PDF Eq. ll3"2T) . 
More generally, using the Lamperti distribution for X[o,t] > 
and probabilistic arguments [33], the PDF of Cta(T', t) 
is approximated by 



Pc TA {T',t>) (z)~[l-P (T,T' 



)] |[1 - P (*', T')] Si l -^dx + [l a (z) + 5 (z)] J 



P (T,V) 



(*) + ¥ 



J381 



where Po(a, 6) is the persistence probability Eq. (I28II . 
Note that to derive Eq. I|38|) we used the fact that 2[o,t] 
and 2T[t',T'] are correlated. In Figs. ElElandElwe plot 
these PDFs of Cta(T' (solid curves) together with 
numerical simulations (diamonds) and find good agree- 
ment between theory and simulation, for the cases where 
these approximations are expected to hold, r > 1/2. In 
the limit t'/T' -» 1 En. ffity simplifies to Eq. QZJ. 



VI. EXPERIMENTAL EVIDENCE 

In this section we analyze experimental data and make 
comparisons with theory. Data was obtained for 100 
CdSe-ZnS nanocrystals at room temperature pflj . We 
first performed data analysis (similar to standard ap- 
proach) based on distribution of on and off times and 
found that a+ = 0.735 ± 0.167 and a_ = 0.770 ± 0.106 
[sol ], for the total duration time T" = T = 3600s (bin size 
10ms, threshold was taken as 0.16max/(£) for each tra- 
jectory). Within error of measurement, a+ « a_ rs 0.75. 
The value of a w 0.75 implies that simple diffusion 
model with a = 0.5 is not valid in this case. An im- 
portant issue is whether the exponents vary from one 
NC to another. In Fig. El (top) we show distribu- 
tion of a obtained from data analysis of power spectra. 
The power spectrum method yields a single exponent 
a psd for each stochastic trajectory (which is in our case 



a+ 



ps 



([)■ Fig- El illustrates that the spread 
of a in the interval < a < 1 is not large. Numerical 
simulation of 100 trajectories switching between 1 and 0, 
with i/>+( T ) = - 0-( T ) an d a = 0.8, and with the same 
number of bins as the experimental trajectories, was per- 
formed and distribution of a values estimated from power 
spectra is also shown in Fig. El (bottom). We observe 
some spread of measured values of a, which is similar to 
experimental behavior. This indicates that experimental 
data is compatible with the assumption that all dots are 
statistically identical (in our sample), in agreement with 

sua- 

We also tested our nonergodic theory and calculated 
distribution of relative on times T + /T, i.e., of the ratios 
of the total time in the state on to the total measure- 
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Figure 13: Histograms of experimental (top) and simulated 
(bottom) fitted values of a for 100 trajectories. Fits are made 
to the power spectral densities of individual trajectories. 



ment time. These relative on times are equivalent to the 
experimental time averaged intensities after their "renor- 
malization" in a way making average intensity in state 
on/off to be 1/0, respectively, in analogy to our model 
stochastic process. Experimental and simulated distribu- 
tions shown in Fig. El are, overall, in good agreement. 
Two important conclusions are derived from these dis- 
tributions of relative on times. First the data clearly 
exhibits ergodicity breaking: distribution of relative on 
times is not delta peaked, instead it is wide in the interval 
between and 1, for different T'. The second important 
conclusion is that for a reasonably chosen threshold (cf. 
Fig- the experimental data is compatible with the 
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assumption 

V>+(t) w ^-(t), 

at least for a wide time window relevant to the experi- 
ments. In other words, not only a + « a_ (ignoring the 
cutoffs) but also A + w A-. This observation cannot be 
obtained directly from the on and off time histograms like 
Fig. because if only power law tails are seen, as in Fig. 
12 these histograms cannot be normalized. To see that 
A + Pi A- note that the distributions of relative on times 
are roughly symmetric with respect to the median value 
of 1/2 (cf. Fig. IHj) . and the ensemble average of relative 
on times is also close to 1/2, while in general the ensemble 
average in our model process is given by A + /(A + +AJ), 
In addition, the variance of the experimental distribu- 
tions for different T" is close to the variance of the Lam- 
perti distribution (1 — a)/4 [33] for a ~ 0.8. There are a 



few comments to make. First, 100 trajectories are insuffi- 
cient to produce accurate histograms, as can be seen from 
the right side of Fig. ideally, these histograms should 
be identical for different X", and given by the Lamperti 
distribution Eq. 13211 . Second, there is an effect due to 
the signal discretization, leading to a natter and wider 
histogram at X" = 36s. Third, there is a certain slow 
narrowing of the experimental histogram as T" increases, 
and the average relative on time slowly decreases. Both 
of these trends are probably due to cutoffs in the power 
law distributions, especially for on times, as can be seen 
in Fig. [2 These trends slightly depend on the choice of 
the threshold separating on and off states. 

As mentioned previously, the groups of Dahan and 
Bawendi 0, li| measure values of a + rj a_ w 0.5 for 
hundreds of quantum dots (see Table QJ, while we report 
on a higher value of a. An important difference between 
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our samples and Dahan/Bawendi groups is that in those 
works the dots are embedded in PMMA, while in our 
case they are not [49] (see also [Tfll|l. 

VII. SUMMARY AND CONCLUSIONS 

Our main points are the following: 

1. Simple three-dimensional diffusion model can be 
used to explain the exponent a = 1/2 observed in many 
experiments. In some cases deviations from a = 1/2 
are observed, and modifications of Onsager theory are 
needed. We cannot exclude other models. 

2. Simple model of diffusion may lead to ergodicity 
breaking. Thus ergodicity breaking in single molecule 
spectroscopy should not be considered exotic or strange. 

3. The time average correlation function is random. 
Ensemble average correlation function exhibits aging. 
Hence data analysis should be made with care. 

4. Our data analysis shows A + fh A_, a + m ct_ (be- 
fore the possible cutoffs) and that the distribution of a 
is narrow. It is important to check the validity of this 
result in other samples of nanocrystals, since so far the 
main focus of experimentalist was on values of a and not 
on the ratio of amplitudes A + /A_. 

How general are our results? From a stochastic point 
of view ergodicity breaking, Levy statistics, anomalous 



diffusion, aging, and fractional calculus, are all related. 
In particular ergodicity breaking is found in other mod- 
els with power law distributions, related to the underly- 
ing stochastic model (the Levy walk). For example the 
CTRW model also exhibits ergodicity breaking j4j, and 
hence a natural conflict with standard Boltzmann statis- 
tics emerges. Since power law distributions are very com- 
mon in natural behavior, we expect that single particle 
ergodicity breaking will be a common theme. Further, 
since we showed that a simple diffusion model can gen- 
erate ergodicity breaking, for the nano-crystals, we ex- 
pect that ergodicity breaking be found in other single 
molecule systems. One simple conclusion is that pre- 
dictions cannot be made, based on ensemble averages. 
In fact the time averages of physical observables remain 
random even in the limit of long measurement time. The 
fact that the time averaged correlation function is a ran- 
dom function, means that some of the experimental pub- 
lished results, on time average correlation functions, are 
not reproducible. 
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